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SUMMARY 


A virial equation of state for nitrogen has been determined by use of 
newly measured speed-of-sound data and existing pressure-density-temperature 
data in a multiproperty fitting technique. The experimental data taken were 
chosen to optimize the equation of state for a pressure range of 0 to 10 atm 
and for a temperature range of 60 to 350 K. Comparisons are made for thermo- 
dynamic properties calculated both from the new equation and from existing 
equations of state. ; 


INTRODUCTION 


Over the past several years, there has been renewed interest in gaseous 
nitrogen in the pressure range of 0 to 10 atm and the temperature range of 60 
to 350 K. This interest has grown out of the development of the cryogenic 
wind-tunnel concept at the Langley Research Center (LaRC) of the National 
Aeronautics and Space Administration (NASA). This new wind-tunnel concept 
uses nitrogen gas at cryogenic temperatures to obtain aerodynamic data at high 
Reynolds numbers and at transonic Mach numbers. The need for an accurate equa- 
tion of state became apparent as real-gas effects on flow simulation and data- 
reduction requirements were studied. Although the well-known equations of 
state (Jacobsen and Beattie-Bridgeman, refs. 1 and 2, respectively) were read- 
ily available, discussions between NASA and the National Bureau of Standards 
(NBS) revealed that very little experimental data were available for gaseous 
nitrogen at the time of their development to substantiate these equations of 
state in the pressure range of interest. Because of the lack of experimental 
data, it was believed that substantial uncertainty could exist in thermodynamic 
quantities derived from either the Jacobsen or the Beattie-Bridgeman equation. 
To eliminate this uncertainty and to provide an equation of state for gaseous 
nitrogen optimized to the more limited pressure and temperature regions used 
with cryogenic tunnels, NBS has taken sound-speed data (ref. 3) for the condi- 
tions of interest and has used a virial equation of state and a multiproperty 
fitting technique (ref. 4) to arrive at an optimized equation of state. 


The authors would like to thank Robert M. Hall and Jerry B. Adcock of 
NASA-LaRC who contributed significantly to the contents of this report. 
SYMBOLS 
Symbols shown below in parentheses are those used in table I. 
a (VEL) sound speed, m/sec 
B (BO) second virial coefficient, liters/mol 


(Bl) first derivative of B 


(B2) second derivative of B 


c (CO) third virial coefficient, (liters/mol) 2 
(G1) first derivative of C 
(C2) second derivative of C 
Cp (CP) specific heat at constant pressure, J/(mol° kK) 
Cy (CV) specific heat at constant volume, J/(mol°K) 
(DIF) percent difference between present values and reference 1 values 
E (E) internal energy, J/mol 
H (H) enthalpy, J/mol 
p (PRES) pressure, atm (1 atm = 101.3 kPa) 
R gas constant, 8.31434 J/(mol°k) 
S (S) entropy, J/(mol°K) 
T ( TEMP) temperature, IPTS-68, K 
V specific volume, liters/mol 
@) (DENS) density, mol/liter 


MULTIPROPERTY FITTING 


The use of sound-speed measurements for calculating thermodynamic proper- 
ties is well established (refs. 5 to 7). The sound speed relates to derivatives 
of the primary thermodynamic quantities Pp, 0, and T ina nonlinear fashion. 
For example, 
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or, equivalently, 


Cp\/ 0 
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It is this nonlinearity which prevents the use of conventional linear least- 
squares fitting techniques (ref. 4). 


The pressure and the sound speed are computed from the density virial 
expansion 


p = REO(l + BO + cp?) (4) 


The virial coefficients B and C are in turn represented by an assumed power 
expansion in temperature as discussed later in the paper. Coefficients from 
which B and Care computed are found by using multiproperty fitting of two 
sets of data. 


Multiproperty fitting is a technique which allows the use of different 
kinds of thermodynamic data simultaneously in a least-squares estimation of the 
adjustable parameters in an equation of state. In the present case this tech- 
nique allows the simultaneous use of p-0-T data and sound-speed data to deter- 
mine a virial equation of state. The technique is especially appropriate for 
this particular combination of data forms to determine the second virial coef- 
ficient B. An additional advantage is obtained in this case because the B in 
equation (4), which is the density virial expansion, and the B in 


pv = RT + Bp (5) 


which is the pressure virial expansion, are interchangeable. Therefore, an 
expression for the sound speed may be derived from equation (5) which contains 
only the variables p, a, and T (i.e., no density). These are the actual 
measured variables. This expression is then combined with equation (4) by 
multiproperty fitting to determine B in the density virial expansion. For 
details of this process, see reference 4. 


EXPERIMENTAL METHOD 


The sound speed was measured with a fixed-path-length cavity, or resonator. 
Electrostatic transducers were used for sound generation and detection. A 
standing wave is produced by exciting the generating transducer with a variable 
frequency. The frequency of the oscillator is adjusted for maximum received 
signal, which occurs at the cavity resonance. At resonance there are an inte- 
gral number N_ of half-wavelengths in the tube of length L, 


The sound speed is computed from 
a = FA (7) 


Generally, frequencies F were measured corresponding to resonances for 

N = 6, 8, 10, 12, 14, or 16. Three to five resonance frequencies were mea- 
sured and the experimental sound speed was taken to be the average of the 
velocities computed from the various frequencies. 


The cavity wall is a tube of Invar 67 mm long and 10 mm in diameter. The 
electrostatic transducers are at the ends of the cavity. The inside wall of 
the tube was ground and polished to allow proper evaluation of the boundary 
effects. The sound speed inside the tube is slower than it would be in an 
unbounded medium because of energy losses to the walls and ends via thermal 
conduction and viscous drag. This difference in sound speed can be computed 
accurately by using thermodynamic properties such as viscosity, thermal conduc- 
tivity, heat capacity, and density of the gas medium (ref. 1). For the present 
experimental method, the size of this correction was generally a few hundredths 
of a percent in sound speed; however, it increased to about 0.5 percent in the 
worst case, which occurred at low temperature and high pressure. 


The active transducer surface was a thin, flexible plastic film (ref. 8), 
coated on one side with aluminum of a thickness of approximately 50 nm (500 IN, 
The conducting side was placed on the ends of the tube and stretched tightly so 
as to produce a flat surface. Immediately on the other side of the coated sur- 
face were placed the electrodes, which were copper discs. The signal strength 
applied to the electrodes was 50 to 100 V at 1 to 30 kHz. The transducers pro- 
duce sound of measurable intensity over this range of frequencies. Over the 
frequency variation of the width of a single resonance the sound intensity is 
essentially constant. 


Temperatures were measured on the International Practical Temperature 
Scale-68 (IPTS-68) with a platinum resistance thermometer calibrated by the 
NBS Heat Section. Uncertainties in the measured temperature were estimated 
to be 5 mK at the lower temperatures and approximately 50 mK at 350 K. Pres- 
sures were measured with a precision quartz spiral Bourdon tube, calibrated 
at the National Engineering Laboratory with an air-dead-weight gage. Pressure 
uncertainties were estimated to be less than 0.013 percent at 15 atm, increas- 
ing to 0.3 percent at 0.33 atm. The sample gas was ultrapure carrier type. 
Analysis of the nitrogen gas by the supplier showed an oxygen concentration of 
1 ppm and total hydrocarbon concentration of less than 0.2 ppm. 


To estimate the uncertainty involved in the overall experimental procedure, 
sound-speed measurements of helium gas were taken in the apparatus and were 


found to be within 0.05 percent of known values. This value of 0.05-percent 
uncertainty is taken to be the uncertainty in the present nitrogen measurements. 


RESULTS 


The density virial coefficients resulting from application of the multi- 
property fit to sound-speed and p-O0-T data are based on temperature expan- 
sions. The form of the second virial coefficient is expressed as 


B= y py rie-j)72 é (GPS O03 ee ee) (8) 
and the form of the third virial coefficient is expressed as 
c=) opr 9/72 (9S 4, 2. Be ww ae US) (9) 


The newly calculated numerical values of B and C are given in the following 
table: 


Virial 
coefficient 


Calculated numerical 
values of B and C 


0.2920642838E-02 
-.2334399080E+00 
. 742508651 8E+01 
-1164300428E+03 
-9253305729E+03 
- 3351 637808E+04 
.1992311 954E+04 
.2189970095E-01 
-.5949349694E-01 
-—.2399796650E+02 

453424441 3E+03 
. 2283997281 E+04 


The values for the thermodynamic properties for nitrogen computed from 
equations (8) and (9) for B and C are given in table I. Values are com- 
puted for isotherms starting at 60 K, increasing at 10-K increments to 110 K 
and at 20-K increments from 110 to 350 K. The second and third density virial 
coefficients with the first and second derivatives of each are shown at the top 
of each page. For most isotherms the property values are computed at pressures 
beginning at 0.2 atm and increasing to 10 atm in 0.2-atm increments. The val- 
ues for two-phase conditions are indicated by the notation "2-PHASE BOUNDARY" 
shown at the right side of the table. Values for the thermodynamic properties 
are given beyond the two-phase boundary as an indication of the properties for 
metastable conditions. 


The tabular values of the thermodynamic properties are arranged for com- 
parison of the values computed from the virial coefficients with the existing 
values from NBS Technical Note 648 (ref. 1). The first column for each case 
is computed from the virial equation. The numbers of the different columns are 
presented as percent of deviation from TN 648 values in the graphs in figures 1 
to 7. However, above approximately 150 K the differences are omitted from the 
graphs because they generally are very close to the values of Technical Note 648 
in that region. 


The largest deviations (approximately 6 percent) were found for c and 
c, for the values at saturation in figure 1 and for Cp for the 100-K iso- 
therm in figure 2. Deviations of about 1 percent were found for internal energy 
(fig. 3) and enthalpy (fig. 4). Deviations of a few tenths of a percent were 
found for entropy (fig. 5), density (fig. 6), and sound speed (fig. 7). The 
largest deviations in property values occurred for isotherms at 90, 100, and 
110 K and at pressures near the two-phase boundary, for which the change in 
sound speed with pressure is largest. 


The sound-speed measurements (ref. 3) from which the virial coefficients 
are computed were taken for isotherms from 80 to 350 K and at pressures from 
0.33 to 15 atm. Figures 8, 9, and 10 show sound-speed deviations at 80, 150, 
and 310 K, respectively, computed with the virial coefficient as the reference. 
These three figures illustrate the accuracy of the three equations - virial, 
Jacobsen, and Beattie-Bridgeman - for predicting the experimental sound speeds, 
and that no one of these equations is clearly superior in terms of fitting the 
sound speed for isotherms over the entire temperature range. However, as 
explained later, the virial equation is expected to be superior in predicting 
specific heat capacities. 


For figures 8 to 10, the largest deviations of predicted sound speeds are 
at 150 K. Sound speeds calculated with the virial equation are 0.1 to 0.2 per- 
cent larger than those calculated from either the Beattie-Bridgeman or the 
Jacobsen equation. For the 80-K isotherm the deviations are in the other 
direction and become larger as the pressure increases, with the largest devi- 
ations occurring for the experimental data. The last three data points are 
believed to be incorrect, since they could not be represented by a virial 
expansion for this isotherm taken separately. However, comparison of virial 
surfaces computed with and without these points shows that the points do not 
affect the derived surface to any degree of significance. At temperatures 
larger than 150 K the agreement with the virial surface was generally as 
seen in figure 10. 


ACCURACY OF THE CALCULATED VALUES OF THERMODYNAMIC PROPERTIES 


The accuracy of calculated values of thermodynamic properties depends on 
a number of different factors. Often an estimation of accuracy is made on the 
basic agreement between the experimental value and the calculated value of the 
property in question. An estimate made in this manner is an estimate of the 
precision of the experimental data and almost always results in an assessment 
of accuracy which is unrealistically small. Systematic errors in experimental 
data are usually greater than the imprecision in those data. 


6 


A more realistic estimate of the accuracy of data is to compare those data 
with similar data from other sources and then make an estimate of the agreement. 
This method is probably conservative, but it is more realistic than equating 
precision to accuracy. In this particular case a realistic estimate of the 
accuracy of a property value may be determined by taking the maximum deviation 
shown for a particular property illustrated in figures 1 to 7. The one excep- 
tion is the specific heat capacity, for which experience has shown that empirical 
equations of state usually predict erroneous heat capacities near the saturated 
vapor boundary (and other boundaries such as the saturated liquid and solid 
boundaries) unless extreme care is taken to prevent the error (refs. 9 to 12). 
This erroneous heat-capacity behavior is due to the high degree of flexibility 
of most equations of state (the Jacobsen equation, for example, has 32 adjust- 
able parameters) and the absence of the data immediately beyond the saturation 
boundary to control the behavior of the equation of state. In the case of the 
saturated vapor boundary, additional strain is placed on the equation of state 
because the equation must evaluate the region of densities between the vapor 
and liquid phases in such a manner as to predict the same pressure given the 
same temperature but different densities (i.e., a van der Waals loop). This 
problem is avoided by choosing the virial equation of state which has a minimum 
of adjustable constants and by using only vapor-phase data (no van der Waals 
loop required) with pressures of 15 atm or less. 


Although multiproperty-fitting techniques used with p-0-T data together 
with either sound-speed data or specific-heat-capacity data can reduce the 
error in specific-heat-capacity values as predicted by empirical equations of 
state near the two-phase transition, neither of these latter two kinds of data 
were available at the time the Jacobsen equation was formulated. 


Based on the preceding arguments, the Cp and cy, from the wide-range 
equations of state near the saturated-vapor boundary are questionable. There- 
fore the maximum uncertainty of the virial and cy is estimated to be 
2 percent, which corresponds to the disagreement shown in figures 2 and 3 
for temperatures for which saturated conditions are not a factor. The average 
uncertainty in Cp and) cy, is..closer to.1. percent. 


SUMMARY OF RESULTS 


For many substances there is a lack of data in the gas region, espe- 
cially near the two-phase boundary. As a consequence, the thermodynamic 
Properties of this region are generally found to be poorly characterized by 
equations of state. To correct this situation for nitrogen, measurements of 
the speed of sound in the gas phase have been combined with p-p-T data by 
multiproperty-fitting techniques to produce a density virial expansion for 
nitrogen. This new virial expansion should provide a more accurate representa- 
tion of the thermodynamic properties (especially heat capacity) in this region. 


The values of the virial coefficients found from this optimization are 
compared in detail to the thermodynamic-property values derived by Jacobsen. 
The virial equation sound speeds are also compared to the new experimental 
data and to those computed from the Beattie-Bridgeman equation and from the 
Jacobsen equation for 80-, 150-, and 310-K isotherms. 


Experimentalists should be encouraged to take more detailed data near 
two-phase boundaries, and especially in the gas region, since the behavior of 
the thermodynamic surfaces near these boundaries is generally the most uncer- 
tain. However, in the case of the gas phase this may well be the most diffi- 


cult region to take data. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 

November 2, 1979 
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Figure 1.- Percent of deviation from TN 648 values at the saturation boundary 
for various thermodynamic properties X. 
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